Stable Three- Axis Nuclear Spin Gyroscope in Diamond 
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We propose a sensitive and stable three- axis gyroscope in diamond. We achieve high sensitivity 
by exploiting the long coherence time of the 14 N nuclear spin associated with the Nitrogen- Vacancy 
center in diamond, and the efficient polarization and measurement of its electronic spin. While 
the gyroscope is based on a simple Ramsey interferometry scheme, we use coherent control of the 
quantum sensor to improve its coherence time as well as its robustne ss against long-time drifts, 
thus achieving a very robust device with a resolution of 0.5mdeg/s/v Hz • mm 3 . In addition, we 
exploit the four axes of derealization of the Nitrogen- Vacancy center to measure not only the rate 
of rotation, but also its direction, thus obtaining a compact three-axis gyroscope. 

PACS numbers: 03.65.Vf, 61.72.jn, 06.30.Gv 



Gyroscopes - sensitive rotation detectors - find wide 
application in everyday life, from GPS and inert ial sens- 
ing, to jerk sensors in hand-held devices and automobiles. 
Conventional gyroscopes are built using micro electro- 
mechanical system (MEMS) technology that allows for 
high sensitivities exceeding 3mdeg/s/\/Hz in a hundreds 
of micron size footprint [lHsj- Despite several advan- 
tages -including extremely low current drives (~ 100/iA) 
and large bandwidths (> 200deg/s)- that have allowed 
MEMS gyroscopes to have ubiquitous usage, they suffer 
from one important drawback: the measurement sensitiv- 
ity drifts after about a few minutes of operation, making 
them unattractive for long-time stable rotation sensing 
required for geodetic applications [4]. The intrinsic rea- 
son for these drifts - the formation of charged asperities 
near the surface of the capacitive transduction mecha- 
nism - is endemic to MEMS, but does not occur in other 
systems, such as atom interferometers [El-Q or nuclear 
spins which are used as gyroscopes. However, to 

achieve sensitivities comparable to MEMS, these systems 
require large sensor volumes (~cm 3 ), long startup and 
averaging times, and large power and space overheads for 
excitation (lasers and magnetic fields) and detection. 

Here we propose to overcome the challenges of these 
two classes of gyroscope by using a solid-state spin gy- 
roscope associated with the nuclear spin of nitrogen- 
vacancy (NV) centers in diamond. Such a sensor (that 
we call nNV-gyro) can achieve sensitivity on the order 
of r\ ~ 0.5mdeg/s/VHz • mm 3 , comparable to MEMS gy- 
roscopes, while offering enhanced stability, a small foot- 
print and fast startup- up time. Thanks to the four pos- 
sible orientations of the NV center, the nNV-gyro can 
operate as a three- axis gyroscope. Furthermore, a combi- 
natorial gyroscope, consisting of a MEMS system coupled 
to the nNV-gyro may allow for the best of both worlds 
- rotational sensing with ultra-high sensitivity and long- 
term stability. 

The nNV-gyro combines the efficient initialization and 
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FIG. 1. A schematic design of the nNV-gyro. A slab of dia- 
mond of 2.5 x 2.5mm 2 x 150/xm is anchored to the device body. 
RF coils and microwave coplanar waveguides are fabricated 
on the diamond for fast control. The NV centers are polarized 
by a green laser (532nm) and state-dependent fluorescence in- 
tensity (637nm) is collected employing a side-collection tech- 
nique 11]. A second set of rf coils rotates with respect to 
the diamond chip frame, for example by being attached to 
one or more rings mounted in a mechanical gimbal gyroscope 
configuration (not shown). 
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measurement offered by the NV electronic spin, thanks 
to its optical polarization and fluorescent readout, with 
the stability and long coherence time of the nuclear spin, 
which is preserved even at high density. 

Similar to NMR gyroscopes, the operating principles 
are based on the detection of the dynamical phase [12j 
that the Nitrogen- 14 nuclear spin-1 acquires when it 
rotates around its symmetry axis. Consider an iso- 
lated spin 1 in a diamond crystal, with Hamiltonian 
Ho = QS% + jNbS z , where b is a small magnetic field, 
7jv& <C Q. The spin is subject to radio-frequency (rf) 
fields in the transverse plane at the frequency Q, on res- 
onance with the <H> ±1 transition. The diamond rotates 
around the spin symmetry axis (z-axis) at a rate Q. In 
a Ramsey sequence (see Fig. [2]) the spin will acquire a 
phase if = (7at6 + Q)t, as the phase of the applied pulses 
changes at a rate ft in the frame co-rotating with the 
diamond. The sensitivity 77 to the rotation rate is shot- 
noise-limited and for unit measurement time is ideally 
given by rj ex 1/ \/tN, where N is the number of nitrogen 



nuclear spins associated with NV centers in the diamond 
chip. The performance of the nNV-gyro is then set by 
the coherence time of the spin, which together with the 
efficiency of its initialization and readout sets the limits 
to the sensitivity, and by its robustness to external pa- 
rameters, such as temperature, stray fields and strains, 
which determines its stability. In the following we ana- 
lyze these effects while presenting details of the nNV-gyro 
operation. 

The nuclear spin is first initialized by polarization 
transfer from the electronic NV spin. The combination 
of a cycling transition on the electronic m s = lev- 
els and a transition through a metastable level for the 
m s = ±1 levels, yields high polarization of the electronic 
spin [H, [lU . The polarization can be transferred to the 
nuclear spin by various methods, including adiabatic pas- 
sage [lj], measurement post-selection Q or exploiting a 
level anti-crossing in the orbital excited state at ~ 500G 
magnetic field @,y,[T3|. It is as well possible to induce po- 
larization transfer at low (or zero field) either in the rotat- 
ing frame J20[ or by exploiting two-photon transitions via 
longitudinal driving [5j. Working at low magnetic field 
simplifies the nNV-gyro operation and avoiding the level 
anti-crossing allows repeated readouts [12| , while still al- 
lowing high fidelity preparation in a time t po i < 2/is (see 
online material for details). A small bias field B ~ 20G 
can be applied to separate the signal from the four differ- 
ent NV classes, associated with the four (111) directions 
in diamond, by shifting them off-resonance. We note that 
the initialization time is much shorter than the startup 
required even for MEMS gyroscopes (a few tens of mil- 
liseconds). 

For ease of operation, we assume that the rf and mi- 
crowave (fiw) pulses used for initialization and readout 
can be delivered by an on-chip circuit, integrated with 
the diamond chip. After preparation, the NV electronic 
spin is left in the |0) state, which does not couple to the 
14 N nuclear spin nor to the spin bath. Then a Ramsey se- 
quence is applied using the off-chip rf driving. A 27r-pulse 
at the center of the sequence, applied with the on-chip rf 
refocuses the effects of stray magnetic fields and provides 
decoupling from the spin bath. We note that the echo 
will be effective even if the on- and off-circuit cannot be 
made phase-coherent. With this control sequence, the 
coherence time of the sensor spin is thus limited by T2 
(and not by the much shorter dephasing time To* ), which 
can be exceptionally long for nuclear spins mm. it 
is then possible to operate at very high density of sen- 
sor spins. Even at a NV density n^v °° 10 18 cm -3 and 
assuming a density of single nitrogen defects (PI cen- 
ters 04261) npi « 10n NV ~ 10 19 cm- 3 , the 14 N T 2 time 
is not appreciably affected by the PI bath. Indeed, while 
the dipole-dipole interaction among PI centers is very 
large at these densities (~ 3MHz), the coupling to the 
nuclear spin is still small (~ 345Hz) due to its low gyro- 
magnetic ratio. This leads to motional narrowing and a 
very slow exponential decay, as confirmed by simulations 
based on a cluster expansion approach |27l-[29|. The 14 N 
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FIG. 2. nNV-gyro control sequence, rf pulses, resonant with 
the quadrupolar transition of the 14 N nuclear spin, are applied 
in a frame rotating at a rate Q with respect to the diamond, 
thus inducing a phase Qt. The nuclear spin is first initial- 
ized by polarization transfer from the NV electronic spin [2?} ■ 
An echo pulse is applied -in the diamond frame- to refocus 
static frequency shifts. Finally, Q is extracted by mapping 
the nuclear spin phase shift onto a population difference of 
the electronic spin, and measuring the corresponding fluores- 
cence intensity. 



coherence time is also affected by the interaction with the 
close-by NV center p], Q , which induces dephasing when 
undergoing relaxation with T\ ~ 2 — 6ms at room tem- 
perature and low field [30j. While in isotopically purified, 
low defect density diamonds at high magnetic field the 
dephasing time T 2 * can be as long 7ms [4J , in the proposed 
conditions of operation we can conservatively estimate 
the coherence time of the nuclear spin to be T2 = 1ms. 
We note that the echo sequence has the added benefit to 
make the measurement insensitive to many other imper- 
fections, such as the temperature variation, strain, back- 
ground stray fields, variation in the quadrupolar interac- 
tion, instability in the applied bias magnetic field, etc. 
Thus this scheme yields a very robust and stable gyro- 
scope. 

After the sensing sequence, the 14 N spin is left in the 
state 
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which can be mapped into a population difference be- 
tween the NV levels, thanks to the A « 2.2MHz hyper- 
fine coupling. The time required to map the state onto 
the NV center is t map = 230ns, which is close to the T 2 * 
time for the NV, thus we expect a reduction in contrast 
due to the NV dephasing. Indeed it is the NV dephas- 
ing time that ultimately limits the allowed spin densities. 
A possible solution would be to perform a spin echo on 
both nuclear and electronic spins to extend the coherence 
time. The nuclear-electronic spin system is described by 
the Hamiltonian 
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AS 2 Z + le bS z + AS Z I Z + QI 2 Z + ln I z , (2) 



where we only retained the longitudinal component of 
the isotropic hyperfine interaction because of the large 
zero-field splitting of the electronic spin, A = 2.87GHz. 
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FIG. 3. nNV-gyro sensitivity in mdeg/s/A/Hz, as a func- 
tion of density. We considered a diamond chip of size 
(2.5 x 2. 5) mm 2 x 150/xm and an interrogation time t = 1ms 
(solid lines) and t = 0.1ms (dashed line). The sensitivity for a 
simple Ramsey scheme (black lines) is limited by the nuclear 
spin T 2 *. Using an echo scheme (red thick line) improves the 
sensitivity, which is now limited by the coherence time of the 
NV electronic spin used to read out the nuclear spin state. 



The readout sequence (Fig. [2]), with pulses on resonance 
to both o ±1 transitions generate the state [31] 

| Vw) = - cos(ftt) |00) + sin(fit) x 

[e <nt (|-l-l) + |+1-1)) + e- <nt (|+l+l> + |-1+1))] /2 

(3) 

Optical readout can then extract the information about 
the rotation Q. We note that the measurement step can 
be repeated to improve the contrast [111, 03 : although 
at low field the effective relaxation time of the nuclear 
spin under optical illumination is short, thus limiting 
the number of repeated readouts when combined 
with high collection efficiency f] m ~ 1 [ll| we can still 
achieve a detection efficiency C ~ 0.25 for n r = 100 
repetitions and a total readout time t ro ~ 150/is [27( ] . 
The higher detection efficiency will also allow to achieve 
a large dynamic-range by exploiting adaptive phase esti- 
mation schemes 0, [33|, Q • The sensitivity can thus be 
estimated taking into consideration all the parameters 
and inefficiencies presented: 
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where we introduce the dead-time td = t ro + t po i, includ- 
ing the initialization and readout time. For a volume 
of V = 1mm 3 and N = n NV V/4, « 2.5 x 10 14 sensor 
spins, the estimated sensitivity for the nNV-gyro is thus 
r] w 0.5(mdeg/s)/\/Hz, better than current MEMS gyro- 
scopes (see Fig. [3j). The stability of the nNV-gyro can 
however be much higher and comparable to atomic gyro- 
scopes. Indeed the echo-based scheme makes the nNV- 
gyro insensitive to long-time drifts due to temperature 
and stray fields. In addition, the NV spin is a sensitive 
probe of these effects, capable of measuring magnetic [Toj| 
and electric [36| fields, as well as frequency [37] and tem- 
perature shifts 138, 39j . The NV spin could then be used 



to monitor such drifts and correct them via a feedback 
mechanism. 

Until now, we assumed that the relative rotation of the 
rf pulses with respect to the diamond was applied along 
the NV axis of symmetry. The nNV-gyro can however 
work as a three-axis gyroscope, extracting information 
about the rotation rate as well as its direction. If the ro- 
tation axis is at an angle (f} with respect to the NV 
axis, the 14 N undergoes a complex evolution since the sec- 
ond rf pulses is not only rotated by an angle ^(i?, (p, fit) 
in the NV x-y plane, but it also has a smaller flip angle 
a(#, ip, fit) < 7r [27|. The state at the end of the Ramsey 
sequence is then 



i(^)-icos(|) cos (?/>)] 



V2 

sin(a/2)cos(V0|0) - 



1+1) 



Q)+icos(f ) cosQ)] . . 
V2 1 i / 



(5) 

The nuclear spins associated with the four different NV 
classes will therefore experience different evolutions; the 
signal from each class can be measured by sequentially 
mapping each class onto the corresponding electronic 
spin via on resonance /iw pulses (a bias field of 10-20G 
is sufficient to lift the frequency degeneracy among the 
4 classes [io|, El| ) . We note that a more efficient scheme 
would take advantage of the repeated readouts and long 
relaxation times of the nuclear spins, to measure the sig- 
nal from 3 NV classes without the need to repeat the 
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FIG. 4. An integrated nNV-MEMS gyroscope, comprising 
a bulk acoustic wave (BAW) single-axis MEMS disk gyro- 
scope Q in a diamond substrate implanted with NV centers. 
The disk acts as a mechanical gyroscope, while the nuclear 
spins contained in it form a spin gyroscope, (a) Schematic of 
the integrated gyroscope operation. The spins implanted in 
the disk are polarized by an on-chip green laser. The elec- 
trodes surrounding the disk are silvered to allow for total in- 
ternal reflection and fluorescence is collected by replacing one 
of the electrodes by an optical waveguide at 638nm. Striplines 
for rf//iw control are fabricated on the disk, (b) Operation 
of the BAW mechanical gyroscope. The BAW is electrostat- 
ically driven in the second elliptic mode through the drive 
(orange) electrodes. A rotation out of the plane causes a de- 
crease in the gap near the sense (grey) electrodes leading to a 
capacitive measurement of the rotation. Combinatoric filter- 
ing with the nNV measurement would lead to noise rejection 
and improved stability. 
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Ramsey interrogation sequence. 

In conclusion, we have proposed a compact solid-state 
device able to measure rotation rates with a resolution 
rj & 0.5mdeg/sec/V / Hz, better than state of the art 
MEMS gyroscopes, while providing improved stability. 
The device takes full advantage of the long coherence 
time of nuclear spins, which is preserved even at very 
high densities, while exploiting its interaction with the 
electronic spin of the NV center for efficient initializa- 
tion and readout. Going beyond our conservative esti- 
mates, the nNV-gyro could achieve the performance of 
inertial-grade gyroscopes with improvements in the co- 
herence time of the nuclear and electronic NV spin, as 



could be obtained with > 50% N to NV conversion effi- 
ciency [iil , with a preferential alignment of the NV sym- 
metry axis along two directions [43|, and with improved 
collection efficiency, exploiting single-shot measurement 
of the NV center at low temperature 0, EH . Alterna- 
tively, one could consider smaller devices -at the micron 
scale- and exploit its long-time stability to improve the 
performance of MEMS gyroscope (see Fig. 2]). 
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Appendix A: Supplementary Material 



1. Polarization scheme for the 14 N nuclear spin 

Polarization transfer between the NV electronic spin 
and the 14 N nuclear spin is complicated by the fact that 
both spins are spin 1. Unlike for spin- 1/2, polarization 
transfer in the rotating frame (under the Hartmann-Hahn 
matching condition) does not lead to perfect polarization. 
Under the high-density conditions here considered, rely- 
ing on selective pulses would require times longer than 
the coherence time of the electronic spin. Solutions to 
this problem have been proposed either working close to 
avoided crossing, where the energy of the nuclear and 
electronic spins are on resonance p]-[3| or by using a 
probabilistic, measurement-based method Both these 
approaches are viable, but have drawbacks. The first 
method prevents the use of repeated readouts while the 
second method is too lengthy. Here we propose to use 
forbidden two-photon transitions to drive the population 
transfer. Driving the NV electronic spin at the A±7&— Q 
transition with a field along its longitudinal (z) axis [El, 
modulates its resonance frequency thus making energy 
exchange with the nuclear spin possible. This is similar 
to two-photon transitions allowed by Floquet theory 
While these transition rates are usually small, the ability 
to drive the NV electronic spin with very high fields [7j, 
makes the polarization time t = tt A ^ b+< ^ short: for a 
Rabi frequency = 500MHz and a field of 20G, the 
time required is only 1.3/is. Given the short dephasing 
time of the NV spin, a two-step process might be re- 
quired, in which the NV is optically re-polarized before 
driving the polarization exchange a second time. 



lattice. The N-to-vacancy axis sets the direction of the 
zero-field splitting axis can then be along any of the four 
tetrahedral (111) crystallographic directions of the dia- 
mond lattice (see Fig|6j). This intrinsic symmetry can be 
exploited to operate the nNV-gyro as a three-axis gyro- 
scope, extracting information about the rotation rate as 
well as its direction. 

While the maximum sensitivity to the rotation rate is 
achieved for the axis of rotation aligned with the sym- 
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FIG. 5. Polarization of the 14 N nuclear spin under longi- 
tudinal drive of the NV electronic transitions, with a Rabi 
frequency Qr = 500MHz and a field of 20G. The simula- 
tion includes dephasing of the electronic spin modeled by 
a Ornstein-Uhlenbeck process yielding a T 2 * time of about 
200ns. To achieve high polarization the process is repeated 
twice by re-polarizing the NV electronic spin (dashed line) via 
optical illumination. 



metry axis, if the rotation Q is about an axis forming 
an angle ip} with respect to the NV axis, the 14 N 
still undergoes a complex evolution that depends on the 
rotation. The second rf pulses is rotated by an angle 
^(i?, ip, Fit) in the NV x-y plane, with 




2. Operation of the nNV-gyro as a three-axis 
gyroscope 

The Nitrogen-vacancy center in diamond consists of 
a substitutional nitrogen adjacent to a vacancy in the 



FIG. 6. Symmetry axes of the NV center. In an ensemble 
of NV centers, the spins will belong to four different families 
(Fi — Fa) depending on their symmetry axis. The rotation 
vector Q will thus induce different dynamics of the four classes 
which can be used to reconstruct the magnitude and direction 
of rotation. 
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The flip angle a($, ip, fit) < tt is also reduced with respect to the nominal angle, 
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16 



2shr 



fit 
T J 



J (cos(2</?) 



1 2 



2 cos(2tf) cos 2 ((/?)) + 3 cos(fit) + 1 



sin 2 (^) sin(2(^) sin 2 



■ cos(#) sin(fit) 



The state at the end of the Ramsey sequence is then 
given by 



V2 

sin(a/2)cos(^)|0) 



icos(-f-) cos(V')] 



l+i) 



e i ^[sin(^)+icos(f ) ( 
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s«0] 



-1), 

(Al) 

where the angles a, ip are different for each family of 
NVs. Measuring the signal from three families thus al- 
lows extracting information about ft. While this scheme 
requires along the transverse direction of each family, us- 
ing a single axis /iw is also possible, although it makes 
the deconvolution algorithm more complicated. 



3. Coherence time of the 14 N nuclear spin. 

We simulated the 14 N nuclear spin evolution under a 
Ramsey and echo sequence. The simulation was per- 
formed using a cluster expansion to model the effects 
of a PI electronic spin bath. The spin bath was modeled 
as an ensemble of electronic spin- 1/2, thus omitting the 
details of the actual bath (such as the strong hyperfine 
coupling of the PI to its nuclear spins, as well as shield- 
ing effects due to nearby carbon- 13 nuclear spins [§| and 
localization effects due to disorder |9[). Still, the simula- 
tions can give rough estimates of the expected T2 times 
and are consistent with the results of a simple model 
based on describing the spin bath as a classically fluc- 
tuating bath described by a Ornstein-Uhlenbeck process. 
In particular, we expect that due to motional narrowing, 
the coherence times are long even at high densities (see 

Pig. ED- 



Improved measurement efficiency by repeated 
readouts 



The detection efficiency is given 
(l + S^)" 1/2 where n 1>0 



by C 
the 



num- 

(no-ni u 

ber of photons collected if the NV spin is in the 
m s = {0, 1} state, respectively. In the repeated readout 



scheme [Tl|, [l2[, the state of the nuclear spin is repet- 
itively mapped onto the electronic spin, which is then 
read out under laser illumination. The measurement 
projects the nuclear spin state into a mixed state, but the 
information about its population difference is preserved, 
under the assumption that the measurement is a good 
QND (quantum non-demolition) measurement. We can 
include the results of these repeated readout by defining 

1-0te 




0.3 0.4 0.5 

time [ms] 



FIG. 7. 14 N coherence decay under a spin-echo sequence 
(black lines) and a Ramsey sequence. We simulated an ensem- 
ble of hundred non-interacting nuclear spin, each subjected to 
a bath of ~ 500 electronic spins, with a density (from bottom 
to top lines) of n = (0.35,1.06,1.76,2.47,3.17) x 10 19 cm" 3 
(The lines for the echo are indistinguishable). 



a new detection efficiency, C Ur = (l + ^(^^2 ) 

which shows an improvement oc -v/tv, where n r is the 
number of measurement. The sensitivity needs of course 
to be further modified to take into account the increased 
measurement time. Provided the time needed for one 
measurement step is smaller than the interrogation time 
(including the initialization time), it becomes advanta- 
geous to use repeated readouts. The maximum number 
of readouts is set by the nuclear spin relaxation under 
optical illumination, driven by non-energy conserving 
flip-flops in the excited states. While at high fields, this 
time is very long, allowing 2000 measurements in 5ms, in 
the present conditions we find that ~ 100 measurements 
would provide a good balance between the improvement 
in C and the reduction due to longer measurement time. 
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